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Using extensive first-principles molecular dynamics calculations, we characterize the superionic 
phase transition and the lattice and electronic structures of the archetypal Type-I superionic con- 
ductor Q-Agl. We find that superionicity is signalled by a phase transition of the silver ions alone. 
In the superionic phase, the first silver shell surrounding an iodine displays a distinct dynamical 
structure that would escape a time-averaged characterization, and we capture this structure in a set 
of ordering rules. The electronic structure of the system demonstrates a unique chemical signature 
of the weakest-bound silver in the first shell, which in turn is most likely to diffuse. Upon melting, 
the silver diffusion decreases, pointing to an unusual entropic contribution to the stability of the 
superionic phase. 
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PACS numbers: 71.15.Pd,66.30.-h,66.30.Dn 



Key advances in energy research have prompted a 
surge of interest in superionic materials, as a crucial 
enabling technology for a variety of nanotechnological 
devices, including sensors, switches, batteries, and fuel 
cells. Of the superionics, Agl and related silver halides 
and sulfides have attracted particular attention because 
of the unusually high levels of ionic conductivity they ex- 
hibit, and as such are findingincreased and varied tech- 
nological implementation 0,12, El- At normal pressure, 
Agl enters its superionic a phase above T c = 420 K, at 
which temperature a phase transition to a body-centered 
cubic structure is accompanied by an increase in the sil- 
ver conductivity of nearly three orders of magnitude, to 
a value of 1.31 tt^ 1 cm -1 Previous molecular dynam- 
ics studies using classical pair potentials have successfully 
reproduced experimental characteristics of the a a nd 8 
phases, as well as the a — > 8 transition [jl lol R, lil 1^. Hof . 
but these are unable to describe the electronic struc- 
ture in a dynamic environment, or to capture the phe- 
nomenology of the melting transition. In this regard, 
first-principles simulations provide unique and unbiased 
predictive power. 

We perform Car-Parrinello molecular dynamics simu- 
lations in the canonical NVT ensemble at temperatures 
ranging from 200 K to 1250 K for a total of 800 ps. All 
simulations were performed with a 54-atom unit cell and 
a = 5.174 A, except the Wannier function calculations, 
which were performed in a 32-atom unit cell [2l| . 

First, we find evidence of a phase transition of the sil- 
ver ions near the experimental T c that is independent of 
the conformation and dynamics of the iodine sublattice 
and signals the transition into the superionic a phase. 
The silvers exhibit a sharp decrease in their diffusion be- 
havior upon cooling below T c , although cubic boundary 
conditions forbid the iodine structural transition to the 
hexagonal wurtzite 8 phase. Results from a series of sim- 
ulations in which we immobilized the iodines in a fixed 
bec configuration provide a secondary, stronger indicator 



of the independence of the silver transition from any io- 
dine dynamics. Fig. ^ displays the associated silver ion 
diffusion coefficients DA g for the fixed-iodine case and for 
a fully mobile lattice. Diffusion coefficients are derived 
from the mean-square displacement (MSD) via the Ein- 
stein relation, D — limt_ > oo (MSD(t)). In both cases, 
the slope of an Arrhenius plot of Dj^ g shows a character- 
istic discontinuity near the experimental transition tem- 
perature. Unexpectedly, this discontinuity is even more 
pronounced for the fixed-iodine case. Immobilizing the 
iodine sublattice does lead to an overall decrease in Da 3 , 
suggesting local lattice fluctuations beneficial for silver 
mobility are frozen out, but the system retains its supe- 
rionic behavior. 

Also of note in Fig. ^ is the unusual decrease of the 
silver ion diffusion coefficient upon melting at about 
850 K, a tendency that has also been observed experi- 
mentally 0, but is not captured by classical poten- 
tials. The thermodynamic implication gives valuable in- 
sight into the stability of the superionic phase. The de- 
creased mobility of the silver ions upon melting signals a 
decrease in the silver entropic contribution, meaning the 
iodines crystallize into the superionic state below T m to 
increase the entropy of the silvers and make them more 
"liquid" . Accordingly, the superionic phase acts as an in- 
termediate between a pure solid and a pure liquid, with 
a high-entropy liquid sublattice flowing through a low- 
energy solid matrix. An entropically driven stabiliza- 
tion of the a phase is also consistent with experimental 
results 0], which determine the entropy difference be- 
tween superionic and liquid phases to be relatively tiny. 

Although it is not possible to separate the total energy 
of the system into respective sublattice contributions, we 
can get a quantitative picture of the energetic fluctua- 
tions associated with each sublattice by instead integrat- 
ing the forces on the ions, which are already decoupled. 
Since in the canonical ensemble, fluctuations a in the to- 
tal energy are related to the specific heat capacity Cy 
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FIG. 1: Arrhenius plot of the diffusion coefficient for the fully 
mobile system (A) and with the iodine sublattice fixed (B). 
The dashed vertical lines indicate the experimental T c (420 K) 
and T m (830 K) at standard pressure. 
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FIG. 2: Specific heat capacities for the silver sublattice, the 
iodine sublattice, and the total system. 

as <t 2 = ksT 2 CV, we can derive a heat capacity- like 
quantity Gy that is sublattice resolved using the ener- 
getic fluctuations derived from the ionic forces. A plot of 
this quantity for the two sublattices, along with the total 
specific heat capacity for the entire system, is shown in 
Fig. [3 Over the full temperature range of the simula- 
tions, we observe the tail end of a divergence in the total 
heat capacity at 400 K and 850 K. These temperatures 
correspond closely to the superionic transition temper- 
ature T c and the melting temperature T m , respectively. 
Examination of the associated ion-resolved C^f 1 curves 
reveals a similar divergence for the silvers at T c but no 
detectable trend for the iodines. As in Fig. ^ we can 
thus link the superionic transition to the silver sublattice 
only. 

We note that an order-disorder transition of the silver 
ions at T c has been reported previously 0, 0] . Our re- 
sults confirm separation of the occupied tetrahedral sites 
into the six inequivalent sublattices of Ref . [15J , the occu- 
pancies of which are plotted as a function of temperature 
in Fig. |3J The simulations reveal an ordering tendency 
for the silvers below T c , characterized by a splitting into 
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FIG. 3: Average silver occupancies of the six inequivalent 
tetrahedral bcc interstitial sublattices of Ref. [lfij . normalized 
against the case for which all ions occupy random tetrahedral 
sites. The inset shows the network of tetrahedral sites sur- 
rounding a single iodine, and the dashed line indicates the 
overall fraction of silvers in tetrahedral sites. 




FIG. 4: (a) The isosurface of silver trajectories at 750 K, 
enclosing a volume associated with the first silver shell sur- 
rounding an iodine, (b) A slice plot of (a), revealing higher 
density (darker) silver occupancy regions. 



higher- and lower-occupancy sublattices. The indepen- 
dent ordering tendency of the silver sublattice upon cool- 
ing is observable despite the inhibition of the a — > (3 tran- 
sition, in agreement with Ref. 0|. However, at these low 
temperatures, even more extensive statistical sampling is 
necessary to allow for comparison of the ordered phase 
with those proposed in the cited studies. 

We have also investigated the most frequented path- 
ways for the silver ions by tracking their positions with 
respect to the iodines and averaging the resulting trajec- 
tories over the course of the simulations. Fig. 0] shows 
isosurface and slice plots for the silver ion density in the 
first shell surrounding an iodine. Our results confirm 
that the highest-density regions lie near the tetrahedral 
sites, with silver ion density smeared toward the octahe- 
dral sites. This agrees well with the experimental results 
of Ref. [Tel. 

From the ionic trajectories, we can identify a set of 
rules governing the instantaneous distribution of silvers. 
The silver-silver pair distribution function gA g -Ag(r) 
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(Fig. |5J) illustrates a zero probability of finding silver ions 
closer together than RA g -A g = 2.4 A, which precludes 
simultaneous occupation of nearest-neighbor tetrahedral 
interstitial sites. Integration of gAg-i(r) over the first 
peak reveals that on average, each iodine has four near- 
est neighbors in the first shell, consistent with geometric 
expectations for tetrahedral interstitial site occupancy. 
In addition, gA g -i(r) indicates the highest probability 
distance for silvers surrounding an iodine is at Raq-i = 
2.6 A, which corresponds to the radial distance to a tetra- 
hedral interstitial site. However, a time-resolved analy- 
sis suggests a more detailed decomposition of the first 
peak of gAg-iif)- We find that most commonly, three of 
the four nearest-neighbor silver ions simultaneously oc- 
cupy a shell that corresponds to the tetrahedral site dis- 
tance (the average value varies from 2.7 to 3.0 depending 
on temperature, with lower temperatures favoring higher 
values). The fourth silver is seen to transition regularly 
between this shell and an equivalent shell for a neighbor- 
ing iodine such that on average, it fills the transition zone 
between the two, as defined by 2.8 < R,Ag-i ^ 4.2 A. The 
transitioning rate is temperature dependent and disap- 
pears rapidly for T < T c . The fourth silver also pos- 
sesses an angular distribution distinct from its three in- 
ner counterparts, as indicated in the inset of Fig. 03 The 
Ag-I-Ag angles for any of the closest three silvers re- 
veal preferences at 65° and 105°, and higher angles are 
surprisingly uncommon. However, the angles introduced 
by the inclusion of the fourth nearest-neighbor shell are 
comparatively diffuse and have significant probabilities 
towards larger values. This suggests that whereas the 
closest three silvers are clustered and correlated in their 
positions, the fourth silver is relatively unconstrained in 
its angular configuration and is affected only marginally 
by the orientations of the remaining three. The rare 
permanent transitions of this unconstrained fourth sil- 
ver between shells represent the driving factor for mass 
diffusion, an observation that would escape experimental 
investigations. 

The above analysis allows for the definition of a set 
of ordering rules that govern the instantaneous distribu- 
tion of silver ions in the first shell surrounding an io- 
dine: (1) four silver ions populate the first shell; (2) no 
two silver ions occupy neighboring tetrahedral interstitial 
sites; (3) on average, three silver ions surround an iodine 
at a radius of Racj-i = 2.6 A, the tetrahedral intersti- 
tial distance; (4) a fourth silver transitions between that 
shell and a second shell at RA g -i 4.2 A, associated 
with a neighboring iodine, and the transition rate be- 
tween the two is temperature dependent and disappears 
below T c ; (5) the angular positions of the three inner sil- 
vers are correlated, with preferred Ag-I-Ag angles of 65° 
and 105°, whereas the fourth (outer) silver is relatively 
unconstrained; and (6) silver ions tend to organize into 
high- and low-occupancy sublattices at temperatures be- 
low T c . 
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FIG. 5: Radial pair distribution functions for (A) Ag-I, (B) 
I— I, and (C) Ag-Ag, along with the integrated Ag-I n(r) 
curve (D), at 750 K. Inset shows the angular distributions of 
the four silver ions in the first shell surrounding an iodine, 
measured with respect to the remaining silvers in the shell 
and resolved according to the three innermost silvers and the 
fourth outer silver. 



We examined the maximally localized Wannier func- 
tions (MLWFs) to obtain a local picture of bonding in 
a dynamic environment |17| . Plotting the time-averaged 
radial distribution of the four iodine Wannier function 
centers (WFCs) about their parent ion reveals an unex- 
pected bimodal separation into short-distance, highly lo- 
calized WFCs and longer-distance, partially delocalized 
WFCs (Fig.|BJ). Further isolating Wannier functions asso- 
ciated with each of the two peaks and plotting the time- 
and statistical-averaged contours for the WFCs around 
the iodines yields the isosurfaces shown in the insets 
of Fig. The long-distance Wannier function centers 
(LWFCs) tend to align along the cubic axes toward the 
octahedral interstitial sites. Their orientations relate to 
the observed smearing of the silver occupancy from the 
electrostatically preferred tetrahedral sites towards the 
octahedral face centers and suggest that these orbitals 
correspond to directional interactions between silvers and 
iodines. On the other hand, the short-distance iodine 
Wannier function centers (SWFCs) exhibit a random an- 
gular distribution, as would be expected in a strictly 
Coulombic picture. Comparing the areas of the two peaks 
in Fig. reveals that of the four WFCs surrounding an 
iodine, 30% on average can be classified as SWFCs, a 
figure which agrees with the likelihood for a first-shell 
silver to be found in the mobile transition zone. This 
value also correlates well with the the experimentally de- 
termined likelihood of a silver to be located outside the 
tetrahedral interstitial site [lH. flif , a quantity confirmed 
by our findings. 

A more detailed picture of the correlation between de- 
localization extent of the iodine MLWFs and the posi- 
tions of nearby silver atoms is offered in Fig. \7\ There is 
clear evidence of a chemical interaction between LWFCs 
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FIG. 6: Histogram of distances between iodine WFCs and 
iodine nuclei. The insets are contours of the iodine WFC 
distribution about their nuclei for WFCs on either side of the 
broken line. 
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FIG. 7: Histograms of interionic bond distance and bond 
angle for the silvers closest to iodine WFCs, resolved into 
(A) long-distance and (B) short-distance Wannier functions. 
The bond angle 9 is defined as the angle between Raq-i and 

RlWFC-I ■ 



and silver atoms within a threshold radius of R,Ag-i ^ 
3.25 A and a solid bond angle of |0| < 20°. This po- 
sitional constraint on the silvers indicates the bonding 
with LWFCs is directional and not merely Coulombic. 
As a quantitative measure of these more complex chem- 
ical interactions, we have also calculated the Born effec- 
tive charges for the ionic sublattices, using the electric- 
enthalpy method of Refs. 19]. We obtain values of 
Z* = ±1.22, larger than the integral values expected 
for a purely ionic crystal, further supporting the conclu- 
sions of our Wannier function analysis. Fig. \7\ shows no 
such directional interaction for SWFCs, indicating that 
corresponding interactions with silvers can be attributed 
to weaker electrostatics. Moreover, examination of the 
radial distribution for first-shell silvers closest to SWFCs 
places them within the transition zone. 

We conclude that silvers contributing to DA g are those 
which are not bound to LWFCs, meaning they are min- 
imally constrained both radially and angularly. Most 
commonly, the outermost of the four first-shell silvers 
fills these criteria, as a lack of stronger directional inter- 
actions with LWFCs increases the average Ag-I bond dis- 
tance. This comparatively unconstrained, mobile fourth 
silver fleetingly occupies the transition zone defined by 



2.8 < Raq-i ^ 4.2 A until it is captured by a neigh- 
boring iodine, leading to mass diffusion. At higher tem- 
peratures, thermal disordering breaks a greater number 
of bonds between silvers and LWFCs, promoting more 
nearby silvers into the transition zone. As such, the over- 
all fraction of occupied tetrahedral sites decreases (Fig.|3J 
and diffusion is enhanced. 

In conclusion, we have shown that the transition to the 
superionic a phase of Agl is signalled by an independent 
phase transition of the silver sublattice alone, character- 
ized by a disordering of the silvers and a sharp increase 
in their diffusivity. Upon melting, the silver diffusion co- 
efficient decreases, pointing to an unusual entropic con- 
tribution to the stabilization of the superionic phase. We 
have also identified diffusion pathways for superionic sil- 
ver ions above T c , and a time-resolved analysis of ion tra- 
jectories has allowed us to define a set of ordering rules 
that govern the instantaneous distribution of silvers in 
the first shell surrounding an iodine. Finally, we have 
found that of the four first-shell silvers, the closest three 
are strongly correlated and restricted in their angular 
distribution, and that they are involved in anisotropic, 
directional bonding to an iodine. The fourth silver is 
bound only weakly and is relatively unconstrained, and 
we have isolated it as the dominant contributor to mass 
diffusion. 
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discussions. 



[1] 
[2] 
[3] 
[4] 

[5] 



[6] 

[7] 
[8] 
[9] 
[10] 



[11] 
[12] 

[13] 
[14] 



40, 
50. 



1337 
1073 



K. Terabe et al., Nature 433, 47 (2005). 

M. Lee et al., Appl. Phys. Lett. 85, 3552 (2004). 

T. Minami et al., Solid State Ion. 86-8, 415 (1996). 

A. Kvist and A. M. Josefson, Z. Naturforsch. 23a, 625 

(1968). 

P. Vashishta and A. Rahman, Phys. Rev. Lett 
(1978); M. Parrinello et al., Phys. Rev. Lett. 
(1983). 

F. Shimojo and M. Kobayashi, J. Phys. Soc. Jpn. 60, 
3725 (1991). 

K. O'Sullivan et al., Phys. Rev. B 43, 13536 (1991). 
F. Zimmer et al., J. Chem. Phys. 112, 6416 (2000). 
J. L. Tallon, Phys. Rev. B 38, 9069 (1988). 
For reviews on the subject, see S. Hull, Rep. Prog. Phys. 
67, 1233 (2004); D. A. Keen, J. Phys. Cond. Matt. 14, 
R819 (2002). 

H. Araki et al., J. Phys. Soc. Jpn. 68, 134 (1999). 

W. Biermann and W. Jost, Z. Phys. Chem. (Frankfurt 

am Main) 17, 139 (1960). 

P. A. Madden et al, Phys. Rev. B 45, 10206 (1992). 
C. Seok and D. W. Oxtoby, Phys. Rev. B 56, 11485 
(1997); Phys. Rev. B 58, 5146 (1998). 



5 



[15] G. Szabo, J. Phys. C 19, 3775 (1986). 
[16] R. J. Cava et al., Solid State Commun. 24, 411 (1977). 
[17] N. Marzari and D. Vanderbilt, Phys. Rev. B 56, 12847 
(1997). 

[18] V. M. Nield et al., Solid State Ion. 66, 247 (1993). 

[19] P. Umari and A. Pasquarello, Phys. Rev. Lett. 89, 157602 

(2002); I. Souza et al., Ibid., 117602. 
[20] S. Baroni et al., URL 



http : //www . quantum- espresso . org/ 

[21] Simulations were performed in a plane- wave basis set us- 
ing the PBE-GGA XC functional, a 4d 10 5s 1 ultrasoft sil- 
ver and 5s 2 5p 5 norm-conserving iodine pseudopotential, 
cutoffs of 22 Ry and 176 Ry for the wavefunctions and 
charge density, and a timestep of 20 au. 



